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Q . Emission of a charged particle propagating in a medium with a curved 

. magnetic field is reconsidered stressing the analogy between this emission 
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mechanism and collective Cherenkov-type plasma emission. It is explained 
how this mechanism differs from conventional Cherenkov, cyclotron or cur- 
vature emission and how it includes, to some extent, the features of each of 



. these mechanisms. Presence of a medium supporting subluminous waves is 



essential for the possibility of wave amplification by particles streaming along 
the curved magnetic field with a finite curvature drift. We suggest an analogy 
between the curvature drift emission and the anomalous cyclotron-Cherenkov 
emission. Treating the emission in cylindrical coordinates in the plane-wave- 



^ , like approximation allows one to compute the single particle emissivity and 

growth rate of the Cherenkov-drift instability. We compare the growth rates 
calculated using the single particle emissivity and using the dielectric tensor 
of one dimensional plasma streaming along the curved field. In calculating 
the single particle emissivity it is essential to know the normal modes of the 
medium and their polarization which can be found from the dielectric tensor 
of the medium. This emission mechanism may be important for the problem 
of pulsar radio emission generation. 
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I. INTRODUCTION 



Studies of relativistic strongly magnetized plasma in astrophysical setting (like pulsar 
magnetosphere) have shown the possibility of a new mechanism of electromagnetic wave 
generation. This new mechanism, which combines features of conventional Cherenkov, cy- 
clotron and curvature radiation, deserves more detailed consideration from the fundamental 
physics point of view. Besides investigating the physical nature of this process, we also 
reconcile different available approaches to this problem. 

In this work we discuss this novel emission mechanism of a charged particle streaming 
with relativistic velocity along curved magnetic field line in a medium. A weak inhomo- 
geneity of the magnetic field results in a curvature drift motion of the particle perpendicular 
to the local plane of the magnetic field line. A gradient drift (proportional to V ■ B) is 
much smaller than the curvature drift and will be neglected. When the motion of the par- 
ticle parallel to the magnetic field is ultrarelativistic the drift motion even in the weakly 
inhomogeneous field can become weakly relativistic resulting in a new type of generation of 
electromagnetic, vacuumlike waves. Presence of three ingredients ( strong but finite mag- 
netic field, inhomogeneity of the field and a medium with the index of refraction larger than 
unity) is essential for the emission. We will call this mechanism Cherenkov-drift emission 
stressing the fact that microphysically it is virtually Cherenkov-type emission process. 

Conventional consideration of the curvature emission ( P|, @]) emphasize the 
analogy between curvature emission and conventional cyclotron emission. To our opinion 
this approach, though formally correct, has limited applicability and misses some important 
physical properties of the emission mechanism. In a separate approach Kazbegi et al. 
considered this process calculating a dielectric tensor of the inhomogeneous magnetized 
medium, thus treating the emission process as a collective effect. They showed that maser 
action is possible only if a medium supports subluminous waves. In this work we show how 
these two approaches can be reconciled and argue that the dielectric tensor approach, which 
treats the Cherenkov-drift emission as a collective process, has a wider applicability. 
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The interplay between cyclotron (or synchrotron) and Cherenkov radiation has been 
a long-standing matter of interest. Schwinger et al. discussed the relation between 
these two seemingly different emission mechanisms. They showed that conventional syn- 
chrotron emission and Cherenkov radiation may be regarded as respectively limiting cases 
of |n — f I ^ 1 and B = of a synergetic (using the terminology of Schwinger et al. [Q) 
cyclotron-Cherenkov radiation. In another work this analogy has been discussed to in- 
clude cyclotron-Cherenkov emission at the anomalous Doppler resonance. An important new 
aspect of our work (as compared with and ) is that we take into account inhomogeneity 
of the medium. 

Physical origin of the emission in the case of Cherenkov-type and synchrotron-type pro- 
cesses is quite different. In the case of Cherenkov-type process the emission may be at- 
tributed to the electromagnetic polarization shock front that develops in a dielectric medium 
due to the passage of a charged particle with speed larger than phase speed of waves in a 
medium. It is virtually a collective emission process. In the case of synchrotron-type pro- 
cess, the emission may be attributed to the Lorentz force acting on a particle in a magnetic 
field. Cherenkov-type emission is impossible in vacuum and in a medium with the refractive 
index smaller than unity. 

Cyclotron emission at the anomalous Doppler effect (cyclotron-Cherenkov emission) is 
an interesting example of Cherenkov-type emission process of a particle in a magnetic field. 
It is impossible in vacuum and requires a superluminous motion of a particle along the 
magnetic field. Thus, the emission at the anomalous Doppler effect is attributed to the 
polarization shock front that a spiraling particle induces in a medium. In our opinion, the 
Curvature-drift emission may be viewed as a Cherenkov-type emission that bears the same 
relation to the conventional curvature emission as the cyclotron emission at the anomalous 
Doppler effect bear to the conventional cyclotron emission. 

In this work we consider a Curvature-drift emission of the particles in the ground gyra- 
tional state. It is possible to obtain the emissivity for particles in excited gyrational state 
by the method of dielectric tensor Then, conventional cyclotron (cyclotron-Cherenkov), 
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Cherenkov and curvature emission mechanisms may be viewed as corresponding limits of 
the Cherenkov-drift mechanism in the cases of homogeneous magnetic field (in a medium), 
medium without magnetic field, and inhomogeneous magnetic field without a medium. 

In Section y we discuss our set up of the problem and how it differs from the previous 
consideration. In Section |ITT| we calculate a single particle emissivity of the Cherenkov-drift 
mechanism and find the growth for kinetic beam instability toward excitation of electromag- 
netic waves at the Cherenkov-drift resonance. Then the results are compared with those 
obtained by the dielectric tensor method. 

II. DESCRIPTION OF A MODEL 

Let us consider concentric coplanar circular magnetic field lines populated with charged 
particles streaming relativistically along the curved magnetic field (Fig. |1|). 
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FIG. 1. Geometry of the considered problem. The magnetic field lines are concentric coplanar 
circular arcs with the radius of curvature much larger than the size of the region. 

The magnetic field is assumed to be very strong, so that the particles are in their ground 
gyrational state (zcroth Landau level). Next we introduce a cylindrical system coordinates 
r, 0, X with the unit vector along the local magnetic field, e^; perpendicular to the osculat- 
ing plane of the magnetic field and Cr = x e^;. We will limit the region under consideration 
to have the size along the r axis br much smaller than the typical radius of curvature of 
the field lines in the region [br <^Rb) and much larger than the typical wave length of the 
considered waves {br ^ A). This physical picture may be considered as an approximation 
to the region of the open field lines of the pulsar magnetosphere. The curved dipole field 
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can be locally approximated by the circular field with the radius of curvature Rb- The 
typical radius of curvature Rb ^ 10^ cm and the wave length of the emitted radiation is 
of the order of meters A ~ 10^ cm. Then for the size of the region satisfying the condition 
lO^cm <^ 5r <^ 10^ cm, it can be considered as homogeneous. 

Next, we assume, that the collective interactions of the particles with waves propagating 
in the system may be described using a dielectric tensor. Then, for a given wave mode, all 
the particles may be separated in the plasma particles, i.e. those particles that interact with 
the wave nonresonantly, and resonant particles. We expect that outgoing electromagnetic 
modes generated by some fluctuating currents at smaller radii will interact with resonant 
particle and will be amplified as they propagate through the interaction region (see Fig. p. 

We wish to find the specific intensities for emission of electromagnetic waves by the 
resonant particles in such a system. This may be done by calculating the work done by the 
extraneous current associated with the resonant particle moving along a given trajectory ( 
^) or by calculating the dielectric tensor of the medium. Below we discuss in details the 
calculations using the former method and compare it with the results obtained using the 
dielectric tensor. 

In a curved magnetic field a particle with relativistic Lorentz factor 7 streaming along 
the field line experiences a drift in the x direction with velocity 



where ub = qB/{mc) is the nonrelativistic gyrofrequency, q is charge of the particle, m is 
its mass and c is the speed of light. 

With the instantaneous 3-velocity v = {O,f0,Md} the radius vector of the particle is 



r(t) = {RB,v^t,Udt}. 



(2) 



The current density associated with the charge is 




(3) 
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We next find the Fourier transform of the current in cylindrical system of coordi- 
nates by expanding the current density in real space (0) in terms of cylindrical waves 
exp {—i{krr + k^x + m(j))) (m - integer) : 

/oo r 
dt / (ir j(r, t) exp (i{ujt — k^r — k^x — m(j))) . (4) 
-oo J 

Using the relation 

/oo 
dtexp (iat) = 2Tc6{a) (5) 
-oo 

and the radius vector we find the Fourier image of the current: 

i{uj,k) = 2nqY exp{-ikrRB}S{uj - k^v^ - k^ua), (6) 

where we introduced /c^ = m/RB- For m 1 we can assume that k^f, is continuous. 

The condition, that the size of the region considered is much larger than the wave length 
places a lower limit on the wave frequency in (^: 

> fi, (7) 

where VL = v^/Rb is the angular frequency of the particle's rotation along the circular 
magnetic field. The condition also vindicates the assumption m ^ 1. 

The expansion in cylindrical coordinates has a limited applicability . Generally, in 
cylindrical coordinates the normal modes will be expressed in terms of Bessel-type functions. 
For the nonresonant modes we can use WKB approximation to the radial dependence of 
normal modes (this is equivalent to the tangent expansion of Bessel functions when argument 
is larger and not close to the order). On the other hand, for the resonant modes the argument 
of Bessel functions is close to the order, so that the WKB approximation (or expansion in 
tangents) is not applicable. In this case we can use Airy function approximation to Bessel 
function, which, in turn, has a wave-like approximation for the interaction of suhluminous 
waves with the particles moving with the speed larger than the speed of light in a medium 
(this corresponds to the Airy function expansion argument is larger than the order). It is 
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shown in that for the large argument expansionof Airy function to apply it is required 
that (i) n — I ^ 1/7^ and (ii) (n — l)m^/'^ > 1 (here n is the refractive index, 7 is the 
Lorentz factor of the resonant particle. Summarizing the above: expansion is valid in 
the two different regimes: (i) for the nonresonant modes, (ii) for the resonant modes when 
n — 1 ^ 1/7^ and (n — l)m?^^ > 1. In what follows we assume that the conditions for the 
expansion (^ are satisfied. 

An important difference of this approach from the one used in and P] is that we 
calculate the transition current in cylindrical coordinates, while in [|l| , and |^ the transition 
current was calculated in Cartesian coordinates. In adopting Cartesian coordinates with 
planar normal modes the interaction length for an individual electron, ^ Rs/'^b, was essen- 
tially coextensive with region over which the waves could possibly interact with any electron. 
This introduces a strict constraint on the particle-wave interaction and precluded a strong 
amplification under all circumstances because the wave would have to grow substantially 
during a very short interaction in a manner that could not be easily quantified. The second 
major difference from and |^ is that we consistently take dispersion into account. 

The major advantage of cylindrical system of coordinates is that the only inhomogeneity 
present in the problem, the weak inhomogeneity of the direction of magnetic field, can be ef- 
fectively eliminated by transforming to the cylindrical system of coordinates. It is the choice 
of cylindrical coordinates that allows one to describe the weakly inhomogeneous system by 
the homogeneous (independent of r ) set of equation, but in the curved coordinates. The use 
of cylindrical waves with the phase dependence of the form exp {i{ujt — kj.r — k^x — mcf))). 
has a limited applicability, but allows for simple estimate of the Cherenkov-drift emissivity 
and growth rate. 

III. CHERENKOV-DRIFT EMISSIVITY AND AMPLIFICATION 

In this section we calculate the growth rate of the Cherenkov-drift instability by finding 
the Cherenkov-drift emissivity for the single particles on the zeroth Landau level streaming 
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along the curved magnetic field. We then compare it with the calculations done using the 
dielectric tensor of a plasma in a curved magnetic field. The results of both calculations, 
obviously, coincide. 

With the known current density, the energy radiated by the particle into a given wave 
mode a is given by ( Q], Eq. (3.18)) 

t/° = 47ri?^(k)|e*"-jK,k)|% (8) 

where e" is the polarization vector of the emitted mode and R% is the ratio of electric to 
total energy in the wave, as defined by @]. 

To find the energy radiated one has first to determine the normal modes of the medium 
and find their polarization. In contract to the vacuum case, where one is free to chose arbi- 
trary combination of plane transverse electromagnetic waves as normal modes, in a medium 
the normal modes must be the eigenvectors of the corresponding dispersion equation. 

To determine the the normal modes and their polarization we must solve the dispersion 
equation 

c2 

Det \kira\ = 0, Aim = {k'^^lm " hkm)^ " Qm, (9) 

where eim is the dielectric tensor of the medium. Dielectric tensor for a one dimensional 
plasma streaming along the strong magnetic field has been calculated in P], [^]. 

Dispersion equation (P) is very complicated - all the nine components of the matrix A;^ 
are nonzero. This is different from the case of a plasma in a homogeneous magnetic field, 
where, due to the cylindrical symmetry around magnetic field, some components of Kim 
could be set to zero without a loss of generality. 

A. Polarization of Waves in Anisatropic Dielectric 

Below we will restrict our consideration to the case of electron-positron plasma with the 
same distribution functions streaming along a superstrong magnetic field. 
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1. Infinitely Strong Magnetic Field 



In the infinitely strong magnetic field the dielectric tensor is 



^ 0^ 



l-K 
1 



(10) 



where 



K 



V4> df^'^ 
to — Qm dpff, 



(11) 



7^ {u — ^Imy 

f{p^) is a one dimensional distribution function. 

The dispersion equation ([TT|) may be factored for the dispersion relations of the t-mode, 
with the electric field perpendicular to the k - B plane, and It-mode, with the electric field 
in the k - B plane (Fig. 0): 

1 



= 1 



nl = 1 



{-n^,0, 



l-K- 



_=_={V(1 -K)il-K- nl)-^, ^(1 - K){1 - K - nl)^} (12) 
^1 — K)^ + Kn']_ '''-L "-L 
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FIG. 2. Normal modes in the considered medium in the hmit of very strong magnetic field. 
The electric field vector of the t-mode is in the plane Cr — ex and the electric vector of the It-mode 
is orthogonal to — k plane. 

These waves are natural analogs of the t- and It-modes in the case of straight mag- 
netic field lines. The t-mode is purely transverse and the It-mode is a mixed transverse- 
longitudinal wave. Generally, in the It-mode the electric field is not perpendicular to the 
wave vector, but for teneous plasma in the high frequency limit if ~ 0, so that the It- wave 
is approximately transverse: 
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2. Finite Magnetic Field 



In the finite magnetic field the dielectric tensor is 



l + 2d 
1-K 
l + 2d 



where for cold plasma d = Up/u]^. Eq. (|^). Dispersion relation is then 



[l-K)ni + {l + 2d){l-K-ni-ni)]{l + 2d-ni-nt-ni) = 



which has solutions 



= (l + 2rf), 



e^*) = —{-n^,0,nr} 



;i + 2d)[l 



1-K 



(14) 



(15) 



(16) 



So that the polarization vectors are the same as in the case of infinitely strong magnetic 
field within factors uj'^/uj\. 



B. Cherenkov-drift emissivity 

With the k vector of the emitted transverse waves given by 

k = {kr-, k(j), kx} , 
we choose the following polarization vectors 



{ kx, 0, kr, } 



(17) 



(18) 



where k± = yk^ + k^. This choice of polarization vectors is a limiting case of very strong 
magnetic field and tenons plasma. 
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We note, that the separation of the normal modes done in m], 0, |^ is related not to 
plane of the real magnetic field, directed along 0, and of the vector k, but to the plane k — ex- 
This is justified only in vacuum or in the homogeneous medium where there a freedom in the 
choice of polarization vectors of the normal modes. In an anizatropic medium the normal 
modes and their polarizations have to be determined from the dispersion equation. 

The eigenvectors ( |18|) are different from those chosen by f^, Q]. The modes chosen 
in those works follow the analogy between the curvature and synchrotron emission. They 
are not the normal modes of the medium. The approach of 0], involving vacuum wave 
polarization and refractive index close but not equal to unity, may be considered as correction 
to the vacuum curvature emissivity due to presence of a medium when effects of the medium 
on wave propagation are small. When effects of the medium on wave propagation cannot 
be considered as small perturbations one has to solve the dispersion relation to find the 
normal modes and their polarization. This stresses once again the importance of a medium 
in what we call the Cherenkov-drift emission. Presence of a medium changes the nature of 
the emitted modes and changes the corresponding emissivities. 

With the polarization vectors (0), the single particle probability of emission (per unit 
volume (ik/(27r)^ may be written as a polarization tensor 



F"^(k,p) 



^ iCA ^ 



C = (19) 

k± 

with a, (3 = t, It. This form of the emissivity may be compared with [Q, Eq. (13.62-13.65), 
and 0. There are two main differences: (i) the approximate expressions in ^ and for 
the single particle emissivity per unit frequency, involving Airy functions, are obtained if 
the transition current is calculated in Cartesian coordinates, while relations ([19|) are exact, 
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(ii) the polarization of the normal modes chosen in and are different from ours. 
The total emissivity, summed over polarizations, follows from (|19]): 



w(k, p) 



kki k I k 



2 



5{uj -k^v^-k^Ud). (20) 



The first term in square brackets corresponds to the emission of the It-mode, second - to the 
emission of the t-mode. 

Next we calculate the growth rate of the Cherenkov drift instability ( [^]): 

T = j dpw;(k,p)^k-^M. (21) 

Growth rate for the It-mode is 



and growth rate for the t-mode is 

r* = ^" / dp^ ( '^^f^ ) ll^^S {u - k^v^ - k^Ud) ■ (23) 



m J \ k±c J dp^ 

Presence of delta function with the Cherenkov- drift resonance condition without the 
gyromagnetic term indicates that this is a Cherenkov-type emission process which requires 
that the medium supports subluminous waves. The necessary condition for the instability is 
also the same as in the conventional Cherenkov instability: the derivative of the distribution 
function must be positive at the resonant frequency and wave vector. In physical terms this 
means that the number of particles with the velocity larger than the phase velocity of the 
waves exceeds the number of particles with the velocity smaller than the phase velocity of 
the waves. This once again stresses the Cherenkov-type nature of the emission. 

It is clear from ( p2D and (|23| ) that the growth rate of the t-wave is proportional to the drift 
velocity and becomes zero in the limit of vanishing drift velocity. As for the It-wave, it can 
be excited in the limit of vanishing drift by the conventional Cherenkov mechanism which 
does not rely on the curvature of the magnetic field lines. We recall, that in the limit of a 
strong magnetic field and oblique propagation It-wave has two branches: one superluminous 
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and one subluminous ( 1TT|]). On the conventional Cherenkov resonance it is possible 
to excite only subluminous waves. We note that the choice of polarization vectors of [|T|, 
and excludes the excitation of the subluminous branch as well (for which electric field is 
not perpendicular to the k vector) thus prohibiting any maser action without taking into 
account drift motion. 

The growth rates (|23|) and (p2|) may be compared with the calculation done using the 
antihermitian part of the dielectric tensor P] . In case of kinetic instability the growth rate 
is given by ( @]) 



(24) 

a;=cij(k) 



;M;^'(ea4/3e/3) 

where e'„^ and e^^ are hermitian and antihermitian parts of the dielectric tensor, co'(k) is the 
frequency of the excited normal mode of the medium, and e is its polarization vector. 

The relevant components of the antihermitian part of the dielectric tensor follow from 

.47rV 2df{p^) 



exx= J dp^u^— b [uj - k^v^ - k^Ud) 



Aii^q^ f 9f(pff)) 
el^= -i J dp^UdV^—^ — 6 {uj - k^v^ - k^Ud) = e'^^, 

Att'^q'^ f dfips) 
^^"^^ ~^~tu~ J ^^'^^'^ dp ~ ^'^^'^ ~ ■ ^^^^ 

Using dl^) and (^5]) we confirm the growth rates (|23|) and (^2]) for ~ c. 

This approach, involving dielectric tensor of the medium may be considered as a more 
general, than the one using the single particle emissivities. The dielectric tensor approach 
takes consistently into account both resonant and nonresonant particles. In calculating the 
single particle emissivities one still has to calculate the dielectric tensor to find the properties 
of the emitted normal modes of the medium. 

Next we estimate the growth rate for the Cherenkov-drift excitation of electromagnetic 
waves in the strongly magnetized electron-positron plasma. In the plasma frame the dis- 
persion relations for the transverse modes in the limit u — k^v^f, ujb for quasi-parallel 
propagation {kr, kx -C A;^) is 
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Lo = kc {1 — 6) , 6 



<1, 



(26) 



(HI); where cu^ = Anq'^np/m is the plasma frequency of electrons or positrons, Up -density 
of plasma, Tp is effective temperature of plasma in units of mc^. 

The resonance condition, given by the delta function in Eq. (^), then reads 

1 



27. 



^ 2kl^ 2\k. c 



0, 



(27) 



cil 



where we used Vres — '^y^ — ^ o : 
resonance is shown in Figs. ^ and ^ 



^). The emission geometry at the Cherenkov- drift 




B 



FIG. 3. Cherenkov-drift emission in the case 5 <^ u^/c^- Drift velocity is perpendicular to 
the plane of the curved field line (B — Rc plane, Rc is a local radius of curvature). The emitted 
electromagnetic waves are polarized along u^. The emission is generated in the cone centered at 
the angle 0"^"^ = Ud/c and with the opening angle (26)^^'^ ^ 9^"^. 
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FIG. 4. Emission geometry of the Cherenkov-drift resonance 



The maximum growth rate for the t-mode is reached when kx/k^ = Ud/c and the max- 
imum growth rate for the It-mode is reached when kr- = 0. We also note, that in the 
excitation of both It- and t-wave it is the x component of the electric field that is growing 
exponentially. 



Estimating (p2D and (^31) using 5-function ( max [ ^ ] ^ cV26/ud and max[ 



ckk I 



c k 



25), we wind the maximum growth rates of the t- and It-modes in the limit 5 ^ l/7res- 



■pc -pi 



UJ 



df{l) \ 



(28) 



where ujp^res is the plasma density of the resonant particles. 

We estimate the growth rates (|28|) for the distribution function of the resonant particles 
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having a Gaussian form: 



V^pt 



1 



exp 



( 



2p? 



) 



(29) 



where ph is the momentum of the bulk motion of the beam and pt is the dispersion of the 
momentum. Assuming in that U(f)b/c ^ 1 we find the growth rates 



where 7;, = pb/ (mc), 7t = pt/ {mc). 

Numerical estimates show, that the growth rate (^) may be large enough to account for 
the high brightness radiation emission generation in pulsars. For the sake of consistency we 
leave the detail investigation of the possible application of this radio emission mechanism to 
pulsar physics for a separate paper. 



In this paper we considered a new Cherenkov-drift emission mechanism that combines 
features of the conventional cyclotron, Cherenkov and curvature emission. We argued, 
that from the microphysical point of view this emission mechanism may be regarded as 
a Cherenkov-type process in inhomogeneous magnetic field. Considering emission process 
in cylindrical coordinates we have obtained the single particle emissivities. We also pointed 
out, that in order to obtain correct expressions for the emissivities it is necessary to use the 
polarization vectors of the normal modes of the medium. Finally, we calculated the growth 
rates of the Cherenkov-drift instability in a strongly magnetized electron-positron plasma. 
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